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ABSTRACT 

We explore the large angular scale temperature anisotropies in the cosmic microwave background due to 
expanding homogeneous local voids at redshift z < 1 . A compensated spherically symmetric homogeneous 
dust-filled void with radius 3 x 10^/z~'Mpc, and density contrast S ^ -0.3 can be observed as a cold spot 
with a temperature anisotropy AT /T ~ -1 x 10"^ surrounded by a slightly hotter ring. We find that a pair of 
these circular cold spots separated by 50° can account both for the planarity of the octopole and the alignment 
between the quadrupole and the octopole in the cosmic microwave background (CMB) anisotropy. The cold 
spot in the Galactic southern hemisphere which is anomalous at the ~ 3(7 level can be explained by such a 
large void at z ^ I. The observed north-south asymmetry in the large-angle CMB power can be attributed to 
the asymmetric distribution of these local voids between the two hemispheres. The statistical significance of 
the low quadrupole is further reduced in this interpretation of the large angular scale CMB anomalies. 

Subject headings: cosmic microwave background - cosmology - large scale structure 
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1. INTRODUCTION 

Despite the success of the cold dark matter (ACDM) model 
in explaining the cosmic microwave background (CMB) 
anisotropy on small angular scales, the possible presence 
of anomalies on large angular scales in the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP; Bennett et al. 2003) data 
has triggered a debate over their origin. 

The reported anomalies are as follows: the anomalously 
low quadrupole observed by COBE/DMR (Smoot et al. 1992) 
and WMAP (Bennett et al. 2003); the octopole planai'ity 
and the alignment between the quadrupole and the octopole 
(Tegmark et al. 2003; de Oliveira-Costa et al. 2004); an 
anomalously cold spot on angular scales ^ 10° (Vielva et 
al. 2004;Cruz et al. 2005, 2006); and an asymmetry in the 
large-angle power between opposite hemispheres (Eriksen et 
al. 2004; Hansen et al. 2004). Evidence for other forms 
of non-Gaussianity on large angular scales has also been re- 
ported (Chiang et al. 2004; Park 2004; Schwarz et al. 2004, 
Larson & Wandelt 2004). It should of course be emphasized 
that the significance of any one of these anomalies is at most 
3 cr. Nevertheless, the accumulation of anomalies suggests 
that there may be hints of new physics to be added to the stan- 
dard cosmological model. The possible implications are suf- 
ficiently important that discussion of theoretical explanations 
is justified. 

Indeed, to explain the origin of the anomalies, various so- 
lutions have been suggested. Luminet et al.(2003) proposed a 
non-trivial spherical topology to explain the low quadrupole, 
and Jaffe et al.(2005) considered a locally anisotropic model 
based on the Bianchi type VII/, universe to explain the 
quadrupole/octopole planarity and alignment. Other papers 
have studied the possibilities that the large-angle CMB is 
affected by local non-linear inhomogeneities (Moffat 2005; 
Tomita 2005a, 2005b; Vale 2005; Cooray & Seto 2005; Rakic 
et al. 2006). None of these explanations is fully satisfactory, 
and it is therefore useful to develop an alternative interpreta- 
tion of the anomalies. 
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Osaka, 577-8502, Japan 

- University of Oxford, Department of Physics, Oxford, OXI 3RH, United 
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In this paper, we explore the possibility that the CMB is 
affected by a small number of local voids at redshifts z < L 
As is well known, the Rees-Sciama effect (Rees & Sciama 
1968) for voids could induce an additional anisotropy on large 
angular scales. Thompson & Vishniac (1987) and Martinez- 
Gonzalez & Sanz (1990) have studied the anisotropy owing 
to a vacuum spherical (Swiss-cheese) void using the thin- 
shell approximation. Panek (1992), Arnau et al. (1993), and 
FuUana et al. (1996) studied voids for a general spherically 
symmetric energy profile based on the Tolman-Bondi solu- 
tion (TBS;Tolman 1934; Bondi 1947). Although the TBS is a 
good description of voids with a general spherically symmet- 
ric density profile, it has a drawback in that it only applies be- 
fore shell crossing. In contrast, the thin-shell approximation 
is suitable for studying voids in the late stages of evolution. 
In section 2, we first derive analytic formulae for dust-filled 
voids in the Einstein-de Sitter background using the thin-shell 
approximation. In section 3, we explore a model of local 
voids that agrees with the observed anomalies on large an- 
gular scales. In section 4, we summarize our results. 

2. DUST-FILLED VOID MODEL 

2.1. Cosmic expansion 

We estimate the temperature anisotropy caused by a homo- 
geneous thin-shell void. For simplicity, we consider a matter- 
dominated Friedmann-Robertson- Walker (FRW) model with 
zero curvature (Einstein-de-Sitter model) as the background 
universe. In what follows, we use units in which the light ve- 
locity c is normalized to 1. In spherical coordinates (r,9,4>), 
the flat background FRW metric can be written as 



ds^ = -dt^ + R\t){dr + r^dn^), 



(1) 



where f is the time and R denotes the scale factor, and dil^ 
is the metric for a unit sphere. Next, we consider a homoge- 
neous spherical dust-filled void with a density contrast S <Q. 
We assume that the size of the void is sufficiently smaller than 
the Hubble radius //"'. Then, the metric of the void centered 
on the origin can be approximately described by the hyper- 
bolic FRW metric as 



ds^ =-dt'^ + R'\t')idr'^+Rlsmh^{r' /Rc)dn^), 



(2) 
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where t' is the time, R' denotes the scale factor, and Rc is 
the comoving curvature radius. Assuming that the curvature 
radius Rc is sufficiently smaller than the void radius r,,, up to 
order 0(_(_H'~^ /Rcf), the scale factor for the hyperboUc FRW 
void can be written as 



R'{t') : 



I) 



2/3 



(f]'//'2)l/3f'2/3 + 3 



lOR}. 



■4/3 



(3) 

In terms of the Hubble parameter contrast Sfj = H' /H-l and 
the density contrast S = p' / p-l, equation (|3} is expressed as 



R' it') ^(1 + 5) 



1/3 



^5(l+(5)i/3l f 



(4) 



where k = H-^/(R'^R^.) = (1 + Sh^ - (1 + 5) is the absolute 
value of the curvature in unit of the background Hubble pa- 
rameter H. Up to order 0(k), equation yields the cosmic 
age for the void in terms of the model parameters Sh and k 
given at a specific epoch as 



1 



1 + Sh 5(1 +5h)- 



2.2. Time solution 



(5) 



In what follows, we calculate the evolution of the internal 
time f ' in terms of the external time t for the expanding void 
(with a peculiar velocity) up to order 0((rv/H~^)^). The co- 
moving radius of the void is expressed as r,,(f) in the exter- 
nal coordinates and as r[,(t') in the internal coordinates. To 
connect the two metrics at the shell, we require the following 
boundary conditions: 

{R{M)f = (R'{t')r[,(t')f ( 1 + ^^^^X^) , (6) 



and 



- dt^ + R\t)dr^ = -dt'^ + R'^(t')dr'^. 



(7) 

where we have assumed that r' <C Rc- Up to order 
0((rv ///"')■'), the curvature term in equation (|6j is negligible. 
Then equation (|6j and Q yield 



dt' = dt[ l-RRr,r\,5H+-R^r^Sl 



(8) 



where a dot means the time derivative d/dt. Let us assume 
that the expansion of the thin shell in the external coordi- 
nates is expressed as r,,(f) cx t^, where /3 is a constant. For 
the Einstein-de Sitter background, the matter dominated uni- 
verse expands as R(t) cx f^/"^. Therefore, equation (|8} can be 
explicitly written as 



dt' = dt\l 



where we define 



77; 



r,m(t) 



Integrating equation (|9j, we find 

2(3 



■t 1 



6/3+1 



18/3+3 



r2 2 



(9) 



(10) 



(11) 



Note that Sh and 77 are functions of t rather than t'. As we see 
in the subsequent analysis, the anisotropy turns out to have a 
leading order of 77^ . Therefore, omitting the curvature term in 
equation (|6j in the process of deriving dt' can be justified. 



2.3. Expansion of voids 

In the Newtonian limit, the peculiar velocity of the void v 
can be written in terms of the gravitational acceleration g and 
the Hubble parameter H and the density parameter f2 as 

2fg 



(12) 



where / is written in terms of the linear growth factor D and 
the scale factor R as 

Equation il3\ can be approximately written as (Peebles 1980) 



0.6 



(14) 



Consider a homogeneous void with an inner radius r_ and an 
outer radius r+. We assume that the mass deficiency inside 
the void -m is equal to the mass of the wall if the contribution 
from the mean density is subtracted. If we assume further that 
the motion of the wall is described by the distance dt, from 
the center of the void where the gravitational mass satisfies 
M(dt,) = m/2, the gravitational acceleration g acting on the 
wall is 



G(m/2) 47r 3 
— , m = —djp, 



(15) 



where p denotes the background density. Plugging g in equa- 
tion il5\ into equation il2\ . il3i . and il4\ . we have (Sakai, 
1995, Sakai et al. 1999), 

v = -]:Hd,:n°-^S. (16) 
6 

Assuming d^, ^ r^, the power of the expansion for a compen- 
sating void is 

/3«-il70"5. (17) 
6 

2.4. Temperature anisotropy 

We denote quantities at the time the photon enters the void 
and those at the time the photon leaves by the subscripts "1" 
and "2", respectively. Primes denote quantities measured by 
a comoving observer in the interior coordinate system and the 
unprimed quantities are measured by a comoving observer in 
the background universe just out of the shell of the void. The 
angles -02' ■02 are defined by reference to figure 1. 

To calculate the energy loss, we apply two local Lorentz 
transformations at each void boundary. The first is to convert 
the photon four- vector momentum in the comoving frame in 
the background universe to the frame in which the shell is at 
rest. The second is to convert it to the frame in the comoving 
frame inside the void. 

The four-vector momentum of the photon that enters the 
void is 

1 

cos 01 

sin0i 




ki=£, 



(18) 



where is the angle between the normal vector of the void 
shell and the spatial three-vector of the momentum. After the 
photon passed the shell, the four- vector is converted to 




(19) 
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Fig. 1 . — Cross section of a spherical void tlirougli tlie plane in wliicli the photon path and the line of sight to the void center lie. Here is the angle between 
the light ray and the normal vector on the spherical void. The subscripts 1 and 2 represent the quantities at the time the photon enters and leaves, respectively. 
Primed variables represent the values inside the void. 



(20) 



2 R[ ' 



(21) 



/ 7i7([l + (^i-^'i)cosV'i-viv'i] 
7i 7i [cos V'l + (vi - v'l) - vi v'l cos ipi ] 
sinV'i 

V 

where vi and v[ are the velocities of the void shell at the time 
t = ti and 7 factors are defined as 71 = 1/(1 - Vj)'/^ and 7} = 
1/(1 -Vj^)'/^. When the photon reaches the far edge of the and 
shell, the four- vector momentum becomes 

1 

cos-!/;2 
sin 


As the photon leaves the shell, the four-vector momentum is 
converted to 

cos 

7272 [ 1 + ( V2 - v'2) COS -02 - V2 V2] 
7272 [cos '1P2 + (V2 - v'2) - V2V2 COS V'2] 



The velocities of the void are 

dry 
Vi = R- 



:E2 



(22) 




dt 
dr' 

dt' 



= f3v(ti) 



■(23) 



(24) 



(25) 



where / = 1,2. From the connection conditions and ([7), 
and equation up to order 0{ry'), vj can be calculated as 



■-?7/(1 + (5h,) 



(26) 



where 77; = 77(f,)and Sm = SH(ti),Ki = (1 + Snif -(I + S(ti)). 
Using the Friedmann equation, the parameters at f = f 1 can be 
written in terms of those at f = f2 as 



1 + S, 



Kl 



1(1' 



^2 { ^2 
U2[W} 



Hf \R 



(27) 



(28) 



From the geometry of the void in the internal comoving frame, 
the relation between the void radius r'^,^ and r'2 is approxi- 
mately given by ^ 

r',[ sin V'l ~ ^",2 sin V'2 • (29) 
The relation between the time f 1 and t2 can be written as 



dt' 
R'(t') 



r^i cos V'l cos V'2- 



(30) 



The energy loss suffered between times f 1 and f2 by a CMB 
photon that does not traverse the void is 



E2\ 



l+Zl 



^1 7 no void 1+^2 



(31) 



where zi and Z2 are redshift parameters corresponding to ti 
and t2, respectively. The ratio of the temperature change for 
photons that traverse the void to that for photons that do not 
traverse the void is 



T ~ V^J void 1+^1 



1. 



(32) 



^ Strictly speaking, equations 1291 and I30i are not exact, because the pho- 
ton path is not generally straight inside the void. However, the curvature 
correction to the photon path is negligible at the order Oiry') (Meszaros & 
Molnar, 1996). 
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Fig. 2. — The temperature anisotropy for a void with r,, = 340/i"'Mpc, 5 = -0.3 and the size/distance parameter vi" = 1 is plotted as a function of the subtending 
angle for various parameters that control the second order effects Left:(e,Q={0,0) (full curve), (0.1,0) (dashed curve), and (0.2,0) (dashed-dotted curve). 
Right:{e, ()=(0,0) (full curve), (0,0.05) (dashed curve), and (0,0.1) (dashed-dotted curve). 



Equations gE, (9),(1 1),(18)-(30), and (32) can be solved re- 
cursively using rj2 as a small parameter. 

After a lengthy calculation, neglecting the terms of order 
O(k^), we find 

At 8 / 

— = -^vl COS i-9pi26H + 4 - «) - ^Shk 

+ 2(Sh(5 -2k) + 7SI+3SI-k) cos^ , (33) 

where 5h = Sniti) and k = (1 +(5//(f2))^-(l + (S(f2)). 

For Sh = -1 and k = 0, we recover the formula for an empty 
Minkowski void (Thompson & Vishniac 1987). Note that the 
formula can also be obtained by considering the hyperbolic 
vacuum limit: Sh = and k = 1 . For a quasi-non-linear void, 
the linear approximation Sh = -(5/3 holds if the background 
universe is matter-dominated (see Appendix). Then equation 
( I33t yields 

AT 



2187 



(5r/^cosV^2 (243/3-((5-42)(5+(5(12+(5)cos(2i/;2) 

(34) 

where S = d{t2). Plugging equation J17> for J7 = 1 into equa- 
tion i34\ . we have a simple formula for a compensating void. 



AT 



5^ril cos 7/;2(3 + 24cos2V'2)- 



2187 



(35) 



Thus, we would observe the dust-filled thin shell void as a 
cold spot surrounded by a slightly hot ring in the sky with an 
anisotropy AT /T « -46^7]^ /El toward the center of the void. 
To generate an anisotropy AT /T = 10"^, for a given density 
contrast 5q at the present, the void radius should be 



-2/3 



3 X 10^ /i-'Mpc, 



(36) 



where //q ' = 3000/!"'Mpc. 

From equation j35> , one can see that the terms that are lin- 
ear in 6 cancel each other if the condition (5// = -i5/3 is im- 
posed. Therefore, we conclude that the anisotropy owing to 
the void is a second-order effect. To model the non-linear ef- 
fects, we introduce two new parameters e and ^ that control 
the second order effect. They are defined as Sh = -Sj'h-eS'^ 
and (5 = -5/6 + (S^ (Morita et al. 1998, Noh& Hwang 2005). 
One can see in the left panel in figure |2] that increasing the 
value of e leads to an enhancement in the redshift of photons 



that enter the neighborhood of the center of the void because 
the increase in the (absolute) curvature k (or equivalently the 
decrease in the density) allows the photons to travel through 
the void in a shorter time. On the other hand, increasing the 
value of ^ leads to an enhancement in the blueshift of photons 
whose trajectories pass near the edge of the void. This is be- 
cause an increase in the expansion rate of the void shell causes 
a larger difference between the gravitational potentials at the 
time that the photon enters and at the time that the photon 
leaves. 

To calculate the angular power spectrum, we need the an- 
gular size 6 of the void. For a void at redshift z centered at an 
angular diameter distance d from the Earth, we have 



Tp(w, 9) = sin '(sin0/w), 



(37) 



where the size/distance parameter is defined as w = R(z)rv/d, 
and 6 is the angular size of the void. From equations i33\ . i37\ 
and non-linear parameters (e,0' can analytically calculate 
the angular power of the anisotropy. 

3. ORIGIN OF LARGE- ANGLE ANOMALIES 

In this section, we construct a model of local voids that can 
explain the observed CMB large-angle anomalies. We assume 
that the voids are "local" at redshifts z < 1 . 

3.1. The quadrupole/octopole alignment and planarity 

Tegmark et al. (2003) have found two new anomalies on 
very large angular scales in the CMB map from the WMAP 
satellite (see also Schwarz et al. 2004; de Oliveira-Costa et al. 
2004). Firstly, the preferred axis of the quadrupole is aligned 
to that of the octopole. Although, the definition of the "pre- 
ferred axis" of these low multipoles is not unique, a simple 
formulation is to think of the CMB temperature fluctuation as 
the wave function 

^(n) = nn) (38) 

and find the axis n around which the angular momentum dis- 
persion 

i(n)= («'|(n-L)2|«') =^m2|fl,„(n)|2 (39) 

m 

is maximized (de Oliveira-Costa et al. 2004). The inter- 
val between the unit vectors lii and 62 corresponding to the 
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Fig. 3. — The octopole concentration parameter f as a function of the 
separation angle B (in degree) between the pair of identical voids. 



"preferred axes" for the quadmpole and octopole is within 
~ 10°(de OHveira-Costa et al. 2004). 

Another anomaly is the unusual planarity of the octopole. 
To measure the octopole planarity, we define the octopole 
concentration parameter T as the percentages of the octopole 
power that can be attributed to the modes with \m\ = 3, 



Tin), 



|fl3-3(n)p+|fl33(n)p 



|a3m(n) 



(40) 



For the TOE map, choosing the preferred direction in which 
the angular dispersion around z-axis is maximal, the octopole 
concentration parameter is T = 94%. The relation of the oc- 
topole planarity and the concentration to the \m\ = 3 modes 
is obvious because the wavenumber for the fluctuation in the 
spherical harmonic 13 -3 or 73 3 as a function of the polar angle 
9 \s,k ^ 1, corresponding to a half- wave on the meridian. 

To explain the quadrupole/octopole alignment and the oc- 
topole planarity, we consider a pair of identical local voids 
with a separation angle between the center 0, = 50° — 60° at 
hi and h2- If the contribution of these cold spots is domi- 
nant over the CMB in the quadrupole and the octopole, then 
the preferred axis for both multipoles is parallel to lii -n2 for 
which the angular momentum dispersion parallel to the ro- 
tation axis is maximum. Therefore, the quadrupole/octopole 
alignment is naturally explained in such a configuration. In 
the coordinate system in which the z-axis is aligned to the 
preferred axis, the equator lies on the equidistant point be- 
tween the two identical voids. Then the wavenumber of the 
fluctuation owing to the two cold spots is roughly ^ ^ 1 for 
the separation angle 6s ^ 60° between the void. As shown in 
figure|3j we find a T = 99.9% concentration for a separation 
angle 6*^ = 53° in the preferred direction independent of the 
size/distance ratio w = Rr^./d. 

We can see in figure|4]that the correlation between the CMB 
quadrupole and the CMB octopole is large only for a cold spot 
centered at {lc,bc) = (-30°, -30°) and a hot spot centered at 
(//,,/?/,) = (-100°, -20°), suggesting that the alignment is due 
to these local fluctuations. Although the hot spot might be 
related to residual Galactic emission, it would be difficult to 
explain the cold spot in this way. For comparison, we have 
simulated the temperature anisotropics for a pair of identi- 
cal voids with radius r,, = 340/i~'Mpc for parameters w = 0.9, 
(e, = (0, 0) and w = 0.47,(e, = (0.3,0). As shown in figure 
13 we can see that the void model can generate a similar cold 
spot as in the observed CMB. The preferred axis is parallel to 
the vector that connects the centers of the voids. 



3.2. The unusually cold spot 

Using a wavelet technique, Vielva et al. (2004) detected 
an unusually cold spot in the Galactic southern hemisphere 
at (Lb) = (-153°, -59°). The significance in the form of 
kurtosis on ^ 4° scales is at the 3(7 (a fluke at the one in 
^ 500) level relative to the simulation of Gaussian fluctu- 
ations. Such a cold spot can be naturally explained by a 
presence of a large void at z ^ 1 in the direction of the line 
of sight to the cold spot. Assuming the Einstein-de-Sitter 
background, for a density contrast S ^ -0.3 at z ^ 1 and co- 
moving radius r,. ^210 Mpc, the temperature anisotropy 

would be AT/T = -46^7]^ /8l 10"^ in the direction of 

the void. The corresponding angular diameter distance at 
z = 1.0 is c/ = 0.29 //q'. Then the angular radius of the void 
would be 6',, ~ 7° and w ^ 0.17. For the flat-A universe 
with ilo = 0.3, the corresponding angular diameter distance 
at z = 1.0 is d = 0.39 //q then the angular radius of the void 
would be 61,, -5°. 

3.3. The north-south asymmetry 

Another anomaly is the so-called north-south asymmetry: 
the systematic difference in the amplitudes of fluctuations on 
opposing hemispheres. Eriksen et al. (2004) first pointed out 
the asymmetry in the power for multipoles 5 < Z < 40 between 
the two hemispheres in which the "north pole" is located at 
(l,b) - (60°, 10°). Subsequent analyses (Hansen et al. 2004) 
confirmed their result. 

In our void model, the asymmetry is naturally explained by 
the asymmetric distribution of voids on the sky. As shown 
in figure |5l the pair of voids at w = 0.9 contributes for mul- 
tipoles / = 2-5 and voids at w < 0.7 contributes for higher 
multipoles / > 5 because the corresponding angular scale is 
I ^ tt/ tan"' w. The number of voids that contribute to the 
CMB large-angle power should be small because the asym- 
metry in the power falls off as 1 /\/N. 

3.4. The low quadrupole 

The origin of the anomalously low quadrupole has been ex- 
tensively discussed since it was first observed by COBE/DMR 
(Smoot et al. 1992). The WMAP team argued that the low 
quadrupole requires a selection at the one in 1 43 level (Ben- 
nett et al. 2003). However, recent studies based on Bayesian 
or frequentist analyses have shown that the statistical signifi- 
cance is actually much less than previously claimed (Tegmark 
et al. 2003, Efstathiou 2004). The key issue is the model- 
ing of the foreground properties. Even though the observed 
power is low, most of the power may be hidden inside the 
Galactic foregrounds. The uncertainty in the amplitude of the 
cosmological background near the Galactic plane will gener- 
ally lower the statistical significance of the observed power. 

In a similar manner, the statistical significance of the low 
quadrupole can be further reduced in our void model. Con- 
sider a pair of voids with identical angular radius 6y = sin" w. 
Then the portion of the sky area masked by the pair of voids 
is approximately / ^ 1 -cos0,,. For w = 0.7, we have 6^ = 44° 
and f ^ 25% which is larger than the portion of the WMAP 
sky masks such as KpO, Kp2 and Kp4 (Bennett et al. 2003). 
Therefore, we expect a further reduction in the statistical sig- 
nificance of the quadrupole if the pair of the voids is located 
sufficiently near to the Earth. 

However, one may object that the amplitude enhancement 
would require an unusual cancellation between the intrinsic 
and the void-induced quadrupoles. This difficulty may be cir- 
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Fig. 4. — The Mollweide projection maps for the cleaned CMB (Tegmark et al. 2003) (left) and those for a pair of identical voids with the size/distance ratio 
w = 0.9 and the non-linear parameters (e.§) = (0,0) (middle) and w = 0.47, (e,^) = (0.3,0) (right). The separation angles of the voids and the void radius are 
assumed to be (9, = 50°, r,, = 340/;"' Mpc, respectively. The gray scale denotes the temperature fluctuations in which the maximum absolute value is set to 50/iK. 
The north pole is aligned to the z-axis for which the angular dispersion .s(ii) of the quadrupole plus the octopole around the z-axis is maximal in the direction 
(l,b) = (-110°, 60°). The coordinate center is located at (l.b) = (-30°, -30°). Top: the quadrupole maps. Middle:the octopole maps. Bottom.the quadrupole plus 
the octopole maps. The boundaries of the void hemispheres are denoted by small light-gray disks. 




2 4 6 8 10 12 14 

I 

Fig. 5 . — Plots of the angular power T; = Q/(/ + l)/27r owing to a pair of identical voids with (e, Q = (0, 0) as a function of / for different locations w = 0.9 
(star), w = 0.7 (box), vv = 0.5 (triangle), w = 0.3 (diamond) for the separation angle 6 = 50° . We assumed S = -0.3. and Vy = 340/i~' Mpc. The distance to the void 
is given by d ~ rv/w at z 0. 



cumvented by the following argument. In the standard ACDM 
model, the adiabatic large-angle anisotropy / < 5 owing to the 
"late" integrated Sachs-Wolfe (ISW) effect (redshift effect for 
a void) partially cancels the "ordinary" Sachs- Wolfe (OSW) 
effect (Hu & Sugiyama 1995). As we have seen in the previ- 
ous section, the photon that traverses a point near the center is 
redshifted. Thus, the presence of the void may provide a fur- 
ther cancellation between the ISW effect and the OSW effect. 



A quantitative estimate is deferred to a later analysis. 

4. SUMMARY AND DISCUSSION 

In this paper, we have explored the temperature anisotropy 
owing to dust-filled homogeneous local voids at z < 1 in 
the flat FRW universe. Local voids with a comoving radius 

~ 3 X 10^h~^ Mpc and a density contrast 6 = -0.3 will give 
a cold spot with anisotropy AT/T -10"^ surrounded by a 
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slightly hotter ring. A pair of such compensated voids at dis- 
tance ~ 4 X 10^/i~'Mpc whose equidistant point in the direc- 
tion (l,b) = (-30°, -30°) can explain the quadrupole/octopole 
alignment and the octopole planarity, and the void at z 1 in 
the direction {l,b) = (-153°, -59°) can be the origin of the un- 
usually cold spot. We also find that the north-south asymme- 
try can be explained by the asymmetric distribution of these 
local voids. 

As we have shown, the temperature anisotropy owing to a 
dust-filled void in the Einstein-de Sitter background is pro- 
portional to S^r^. Assuming the CDM power spectrum for 
the linear perturbation, the expected temperature anisotropy 
owing to the void is maximum at the scale r,, ^ 100/;~'Mpc, 
because S cx for the linear regime whereas 5 cx '''^ for the 
non-linear regime. Thus the contribution from the single non- 
empty voids at the scale r,, ^ 100/i~'Mpc is more important 
than the empty voids at smaller scales. 

For simplicity, we have not considered the effect of the cos- 
mological constant. However, we expect that the change in 
the deceleration parameter will not dramatically alter the or- 
der of magnitude of the anisotropy although the presence of 
the linear ISW effect ~ 0{Srl) may somehow affect the fac- 
tor. The amplitude owing to the second-order effect would be 
^ 0(5^rl). Even if the density contrast 6 is small, the corre- 
lation between the linear ISW term and the second order term 
can enhance the observed signal (Tomita 2005a,2005b). The 
extension of our analysis to a flat-A FRW background will 
be presented in the forthcoming paper (K.T. Inoue & J. Silk 
2006, in preparation). 

We did not go into details about the non-linear evolution 
of the voids. Our thin-shell approximation should provide a 
good description of voids that are sufficiently smaller than the 
Hubble radius in the late stages of evolution. However, in 
the middle stages, the density fluctuation and the peculiar ve- 
locity of the shell may depend on the initial configuration of 
the density and the velocity perturbation. As we have seen, 
the temperature anisotropy owing to a quasi-non-linear void 
is a second-order effect (see also Seljak, 1996; Cooray 2002; 
Tomita 2005a, 2005b), which may significantly enhance the 
anisotropy. Therefore, we need a more elaborate calculation 
to make a more precise prediction about the anisotropy at- 



tributed to the local voids. 

We have not addressed the issue of void origin. One may 
find discussions of this in previous treatments of the impact 
of voids on small angular scale CMB observations (Griffiths 
et al. 2003). The voids postulated here represent the large- 
radius tail of the void distribution discussed in these earlier 
studies. They represent a non-Gaussian fluctuation distribu- 
tion that is linearly superimposed on the usual scale-invariant 
spectrum of Gaussian adiabatic density fluctuations. Mathis 
et al. (2004) discussed the implications of a distribution of 
primordial voids for galaxy and cluster formation. There are 
similar implications for the large voids that we have postu- 
lated if these are part of a statistical ensemble of primordial 
voids. In rare regions, associated with the thin shell surround- 
ing the voids, cluster and galaxy formation will be enhanced. 
In principle this will be detectable via massive clusters and 
galaxies at redshifts significantly larger than expected in the 
purely Gaussian model. Examples include detection of mas- 
sive galaxies and galaxy clusters at high z and massive clusters 
with high concentrations of dark matter. All of these cases, if 
confirmed, would represent potential difficulties for the stan- 
dard model of primordial Gaussian fluctuations. 

Another implication unique to the very large ^ 300 Mpc 
voids postulated here concerns increased dispersion in the 
Hubble constant as measured both in different directions and 
at different redshifts (Tomita 2001). The dispersion will be 
positively skewed since voids tend to generate a larger Hubble 
constant. Moreover, the global value of the Hubble constant 
will be lowered. The effect could be as large as ^ 10%. 

So far, the large scale voids we have assumed have not been 
detected. In the future, the ongoing projects such as the 6dF 
galaxy survey (Jones et al. 2004) may detect our postulated 
voids (underdensity region). 
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APPENDIX 

PERTURBATION OF HUBBLE PARAMETER 

Let us write the Friedmann equation at a point x in a perturbed universe in terms of the scale factor a, the Hubble parameter H, 
the gravitational constant G, and the curvature K as 



H\x,t)=—pix,t) 



K(x,t) 



3 a\t) 



(Al) 



where the Hubble parameter H is defined on a comoving slice. >From the perturbed equation iAH in the Fourier space, for the 
flat FRW background without A, we have 

2 

(A2) 



where the curvature perturbation TZ is defined as 



(A3) 



From the Poisson equation, the Newtonian gravitational potential $k is written in terms of the density contrast (5k as 



(A4) 
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For a constant equation of state w ^-1, the continuity equation gives a relation between the gravitational potential <i>k and the 
curvature perturbation TZk as (Liddle & Lyth 2000) 



3 + 3w^ 



From equations (IA2> . JA4> . and (IA5> . we have 



H 



5 + 3w 
4 



(A5) 



(A6) 



3(1 +w) 

Thus, for the matter-dominated flat FRW universe w = 0, the Hubble parameter contrast is written as Sh =-S /3 in the linear order. 
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